Coulomb Gap: How a Metal Film Becomes an Insulator 



V.Yu. Butkofcl, J.F. DiTusa, and P.W. Adams 
Department of Physics and Astronomy 
Louisiana State University 
Baton Rouge, Louisiana, 70806 
(February 1, 2008) 

f~* > ) ■ Electron tunneling measurements of the density of states (DOS) in ultra-thin Be films reveal 

that a correlation gap mediates their insulating behavior. In films with sheet resistance R < 50000 
the correlation singularity appears as the usual perturbative ln(V) zero bias anomaly (ZBA) in the 
DOS. As R is increased further, however, the ZBA grows and begins to dominate the DOS spectrum. 
This evolution continues until a non-perturbative \V\ Efros-Shklovskii Coulomb gap spectrum finally 
emerges in the highest R films. Transport measurements of films which display this gap are well 
described by a universal variable range hopping law R(T) = ( h/2e 2 )exp(T /T) 1/2 . 
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PACS numbers: 71.30.+h, 72.15.Rn, 73.40.Gk 

It has been known for some time now that 2D is the lower critical dimension for disordered transport and that even 
a non-interacting 2D electron gas will be localized in the presence of arbitrarily small disorder in the thermodynamic 



, limit [yj2|. When these systems are probed at a finite length scale, by a magnetic field for instance, logarithmic 
corrections to the Drude conductivity are seen but one does not expect a true metal- insulator transition in 2D or 

Seven a metallic phase in the conventional sense 0]. What is somewhat less well understood are the ramifications 
Coulomb interactions and their attendant correlations, particularly in 2D systems with moderate to strong disorder. 
This issue in particular has now come under considerable scrutiny with the recent discovery of an apparent metal- 
insulator transition in the dilute 2D electron gas of Si metal-oxide-semiconductor field effect transistors (MOSFET's) 
|5],Q. This somewhat surprising discovery runs so counter to conventional wisdom that there has been speculation 
that perhaps e — e interaction effects are stabilizing an anomalous 2D metallic phase . Clearly, a systematic study 
of the DOS spectrum and corresponding transport characteristics of an increasingly disordered 2D electron system is 
needed to help illuminate the crucial interplay between disorder and correlations as the system is brought from the 
weakly to the strongly localized regime. In the present Letter we present such a study using ultra-thin Be films with 
■ low temperature sheet resistances ranging from R = 50017 to 2.6MQ. 

The theory of interaction effects in disordered electronic systems has for the most part been developed in two 
extreme limits. In the weak disorder/interaction 2D limit it is known that the primary effect of e — e interactions is 
to produce a logarithmic suppression of the density of states (DOS) at the Fermi energy B], SN ~ — ln(V). This is 
commonly known as the zero bias anomaly (ZBA) and has been well established in a number of different systems via 
tunneling measurements of the DOS This depletion of the DOS is perturbative and results in a weakly metallic 

ln{T) transport conductivity In the opposite limit, i.e. strongly insulating regime, Efros and Shklovskii M]ll| 
have shown that the Coulombic interactions can produce a non-perturbative gap in the DOS, which is commonly 
known as the Coulomb gap. The 3D Efros-Shklovskii Coulomb gap, which has a quadratic energy dependence, has 
only recently been observed Jl2],[l3| . Interestingly, the 2D Coulomb gap is expected to be linear in energy , 

3 i N(eV) = leVl , (1) 



. £h ' where k is the relative dielectric constant, e Q is the permittivity of free space, and a is a constant of order unity. Despite 
Eq.(l) having been in the literature for more than 20 years now, there has been no direct spectroscopic verification 
of its peculiar linear energy dependence due to the intrinsic technical difficulties associated with measuring the DOS 
in thin insulating films pi} ]. Nevertheless a direct measurement of this gap is important in that it is believed to play 
a crucial role in the electron transport |l5[ ] and thermoelectric |l6| properties of highly disordered 2D metals and 
semiconductors. For instance, it is well known that insulating films typically obey a modified variable range hopping 
law of the form 

R(T) = R e X p(T /T)», (2) 

where R is the film sheet resistance, and R is a constant. In the case of a flat DOS near the Fermi energy, the film 
simply obeys Mott's variable range hopping law with v = 1/3 |T^]. If there is a simple gap in the DOS then T Q is 
the gap energy for fixed range hopping and v = 1. Finally, if the DOS spectrum is given by Eq.(l) then one expects 
v = 1/2, R to be of the order of the quantum resistance Rq = h/e 2 Q and 
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T = 2.8e 2 /(fcs47re O, 



(3) 



where ks is the Boltzman constant and £ is the localization length pi . Thus when the v = 1/2 hopping form is 
observed, the DOS spectrum is simply assumed to be that of Eq.(l). In this Letter we examine this assumption 
by presenting the first systematic electron tunneling study of the DOS in uniformly disordered metal films whose 
transport properties range from that of weakly metallic to strongly insulating. We show that with increasing film 
resistance the DOS evolves from an essentially flat spectrum with a perturbative ln(V) ZBA to that of a Coulomb gap 
given by Eq.(l). We also show that the emergence of the Coulomb gap coincides with the emergence of the v = 1/2 
hopping behavior of Eqs.(2) and (3). 

The Be films used in the present study ranged in thickness from 1.5 - 2.0nm with corresponding sheet resistances 
R = 500£7 — 3Mtt at T = bOmK. They were deposited by thermally evaporating 99.5°/« pure beryllium powder onto fire 
polished glass substrates held at 8AK. The evaporations were made in a 4xl0 -7 Torr vacuum at a rate ~ 0.30nm/s. 
The film area was 1.5mm x 4.5mm. Scanning force micrographs of the films' exposed oxide surface did not reveal 
any salient morphological features down to the 0.7nm resolution of the instrument. In fact, the films seemed to be as 
"smooth" as the fire-polished glass on which they were deposited. In addition, a transmission electron microstructural 
analysis of 15nm thick Be films deposited on cleaved NaCl crystals at 84K revealed that the films were composed of an 
ultra-fine base structure that was interspersed with 5 — I5nm Be nanocrystallites. Electron diffraction measurements 
showed no diffraction from the metallic base structure suggesting that it was amorphous. Similarly the oxide (BcO) 
produced a broad, continuous diffraction ring indicating its grain size was < lnm |l9f| . 

The tunnel junctions were formed by exposing the films to atmosphere for 0.1 - 3 hours in order to form a native 
oxide, then a 20nm thick Ag counter-electrode was deposited directly on top of the film with the oxide serving as 
the tunnel barrier. The junction area was 0.7mm x 0.7mm. This technique produced tunnel junction resistances of 
order Rtj ~ Ikfl to 1000/cf2 depending upon the exposure time and other factors. We were always careful to ensure 
that Rtj >> Rfu m at low temperatures. Our lowest resistance Be films superconducted ]Ts| which allowed us to test 
the integrity of the junctions by measuring the dc I-V characteristics at temperatures well below the superconducting 
transition temperature T c . The sub-gap impedance of a "good" junction was always greater than 10 8 £1 |2C|] . 

Shown in the inset of Fig.l is the normal state conductance of the 2600S1 film as a function of ln(T). This sample 
had a superconducting transition temperature T c = 0.33K which was suppressed by the application of a magnetic 
field oriented along the film surface. (T c was a monotonically decreasing function of R.) This sample had the expected 
ln(T) weakly insulating behavior with some rounding below lOOmK. In stark contrast, the main body of Fig. 1 shows 
the activated-like behavior of the 2.6MH film. The solid line is a least squares fit to the data. The linearity of the 
data in Fig.l shows unequivocally that the hopping exponent v — 1/2. Furthermore the slope and intercept of the fit 
determine the parameters T Q ~ 1.6K and R ~ Rq/2 in Eq.(2). It appears that R a is universal in that it is of order 
the quantum resistance. 

In principle we can also use Eq.(3) along with the measured value of T Q to calculate a localization length £. However 
the relative dielectric constant k in Eq.(3) is unknown for a highly disordered metal film. Alternatively, if one instead 
assumes a reasonable value of £ ~ lnm and takes T a ~ IK then Eq.(3) predicts n ~ 10 4 . This value seems reasonable 
in that it lies between the metallic and insulating limits of k ~ oo and k ~ 10 respectively. This is not an issue in 3D 
semiconducting systems where the dielectric constant is well defined in the insulating phase. In any case, it is evident 
that Fig.l is consistent with the existence of the Coulomb gap described by Eq.(l). Accordingly, if the Coulomb gap 
is indeed measurable in 2D, then it will be seen in a sample such as that of Fig.l. 

Electron tunneling provides one of the most direct experimental measures of the DOS spectrum of an electronic 
system. The usual tunnel junction geometry consists an electronically "inert" metallic counter-electrode which is 
separated from the film by a highly insulating barrier, usually a few angstroms of oxide. When a potential is placed 
across the barrier electrons tunnel from the counter-electrode and into the film. At low temperatures the tunnel 
current is proportional to both the counter-electrode's and the film's DOS For the technique to be useful in 

highly resistive films, however, care must be taken to insure that the impedance of the tunnel junction is much higher 
than that of the film. Otherwise, some of the measured voltage drop will actually occur in the film itself, since the 
tunnel current must be drained off through the film. Additionally, there is a legitimate concern that screening effects 
associated with the close proximity of the highly conducting counter-electrode might "wash out" any Coulombic 
features in the DOS. Fortunately, this has been shown not to be the case in recent measurements of the Coulomb gap 
in3DSi:B@. 

In order to demonstrate the resolution of our tunnel junctions and the 2D nature of the Be films, we have plotted 
in Fig. 2 the tunneling conductance of a R — 53051 superconducting film, T c = 0.55K, in a parallel magnetic field 
just below its parallel critical field H c ii = 1.2T. Note that we can resolve two peaks on either side of V = 0. These 
peaks result from the Zeeman splitting of the usual BCS DOS peak and indicate that the superconductivity was 
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spin-paramagnetically limited [ |22| . The inset of Fig. 2 shows the normal state tunneling conductance of the same 
film in a supercritical parallel field. The normal state displays the expected ln{V) DOS anomaly which is consistent 
with the the ln{T) transport behavior shown in the inset of Fig.l. 

In Fig. 3 we show the evolution of the DOS with increasing disorder by plotting tunneling spectra of Be films with 
T = 50mK sheet resistances of R= 530ft, 2600ft, 16000ft, and 2.6Mft. The 530ft and 2600ft samples superconducted 
with T c = 0.55-fT and T c = 0.33K respectively. For those samples a parallel magnetic field was applied in order to 
suppress the superconducting state and a standard ac lock-in technique was used to measure tunneling conductance 
G directly. In the higher resistance films, the conductance was obtained by numerically differentiating dc I-V curves. 
Note that in each data set there is a significant ZBA. In the 530ft and 2600ft films the ZBA was of the weak disorder 
form 5N/N ~ — ln(V). However, as is clearly evident in Fig. 3, the ZBA was no longer perturbative in films with 
R > 10 4 ft. An exponential growth in the ZBA is shown in Fig. 4 where we have plotted G(0)/G(l5mV) as a function 
of R. By making a linear fit to the data we find that G(0) ~ exp(—R/R c ) where R c = 6000ft » Rq/4. 

The most interesting attribute of the data in Fig. 3 is that the conductance traces become more linear at the highest 
resistances studied. The energy dependence of the 16000ft curve represents the DOS spectrum of a film somewhere 
between the weak and strong disorder limits and is correspondingly neither ln(V) nor \V\ in form. In fact, this 
spectrum is best described by a V 1 / 2 dependence. In contrast, the 2.6Mft spectrum is quite linear and symmetric. 
For a consistency check, we also measured the dc I-V characteristic of the films over the same current and temperature 
ranges that were used in the tunneling measurements. The nonlinearity in these transport I-V's were negligible in 
comparison to the tunneling I-V's indicating that the traces in Fig. 3 are representative of the tunneling behavior. The 
solid line through the 2.6Mft data is a least squares fit to the form G{V) — (3\V\ where f3 is an adjustable parameter. 
We believe that the linear behavior of this data represents the Efros-Shklovskii gap as given by Eq.(l). Unfortunately, 
we were unable to determine an absolute normalization of the tunneling conductance. The spectrum remained linear 
up to our maximum biases of ~ 100mV\ Consequently we cannot make a quantitative comparison between Eq.(l) 
and the tunneling data. Nevertheless, the observation of a linear spectrum in films that show the expected hopping 
transport behavior is compelling. 

In conclusion, we have made the first direct spectroscopic measurements of the Efros-Shklovskii Coulomb gap in 
a 2D system. We find the the gap emerges from an exponential growth in the ZBA as the film sheet resistance is 
increased. The details of the activated transport behavior of films showing the gap are in excellent agreement with 
theory and have a universal prefactor of the order of the quantum resistance. In principle, spectroscopic studies of the 
Coulomb gap in films deposited in situ would allow for an absolute calibration of the tunneling conductance since a 
variety of film thickness could be investigated using a single tunnel junction. Such measurements could then be used 
in conjunction with transport measurements to extract the two primary microscopic parameters of the theory k and 
£, via Eqs.(l) and (3). 
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FIG. 1. Semi-log plot of the resistance of the 2.6MQ Be film as a function of T -1 ^ 2 . The solid line is a linear fit to the data 
for which we get T = 1.6K and R a w h/(2e 2 ), see Eq.(2). Inset: Conductance of the 260057 film as a function of ln(T). The 
solid line is a guide to the eye. 

FIG. 2. Tunnel conductance in the superconducting state of the 530fi film in a parallel magnetic field H\\ = LIT. The 
critical parallel field was H c \\ = 1.2T. The tunneling spectrum demonstrates that we can resolve the Zeeman splitting of the 
usual BCS DOS. Inset: Normal state tunneling conductance of the same film as a function of ln(V) at H\\ = 2.5T. 



FIG. 3. Tunnel conductances normalized to G(15mV) for Be films with T = 50mK resistances of R — 530f2, 2600f2, 16000J1 
and 2.6MQ, (top to bottom). The solid lines are a best fit to the form G(V) = (3\V\, where /3 is an adjustable parameter. The 
2.6MQ film had a R T j = 0.6Mfi junction. To insure that R fam « Rtj we took the 2.6MQ data at 700mK, see Fig. 1. The 
other curves were measured at 50mK. 

FIG. 4. Semi-log plot of the normalized zero bias tunneling conductance as a function of 7?. The solid line is linear fit to the 
data. 
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